Concavity of the quantum body for any given dimension 
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Let us consider the set of all joint probabilities generated by local binary measurements on two separated 
quantum systems of a given local dimension d. We address the question of whether the shape of this quantum 
body is convex or not. We construct a point in the space of joint probabilities, which is on the convex hull of 
the local polytope, but still cannot be attained by measuring d-dimensional quantum systems, if the number of 
measurement settings is large enough. From this it follows that this body is not convex. We also show that for 
finite d the quantum body with POVM allowed may contain points that can not be achieved with only projective 
measurements. 
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Quantum correlations are points in the space of correlations 
which are achievable in quantum phyics by performing lo- 
cal measurements on separate quantum systems. In contrast, 
classical correlations can be achieved by local strategies us- 
ing shared randomness. For a given number of measurement 
inputs and outputs the set of classical correlations forms a con- 
vex polytope 111]. However, we have learned from the theorem 
of John Bell that there exist quantum correlations that lie out- 
side this polytope ||2tl. Thus the set of quantum correlations 
(which we refer to as the quantum body) is strictly larger than 
the set of classical correlations. 

Let us first consider the quantum body consisting of two 
spacelike separated parties, each having a choice of perform- 
ing two measurements with two outputs. If we keep only joint 
correlations in the set (excluding marginal terms), we obtain 
the simplest nontrivial quantum body. The boundary of this 
set has already been described by Tsirelson [3], notably, de- 
riving the maximal quantum violation for the Clauser-Horne- 
Shimony-Holt (CHSH) inequahty |4]. Subsequent works lH] 
characterized the boundary of this quantum domain in differ- 
ent but essentially equivalent ways. Recently, the structure of 
this body has been the subject of analytical study in Ref. ^ 
deriving quadratic inequalities. However, the fact that in this 
quantum body marginal terms are not included, the inequali- 
ties derived can give only partial information on the full prob- 
ability distribution, i.e., on the shape of the whole quantum 
body for two parties with two inputs and two outputs. 

Beyond this scenario, Navascues et al. Q devised a so- 
phisticated method based on a hierarchy of semidefinite re- 
laxations. This is completely general, in that it can be ap- 
plied to any number of parties, performing measurements with 
any number of inputs and outputs. However, the method in 
its present form works efficiently for the case when no con- 
straints are imposed on the dimension of the system. In ab- 
sence of such a powerful program, the shape of the quantum 
body for a fixed dimension is not well understood. It is not 
even known whether it is convex or concave. Without the re- 
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striction for the dimensionality of the quantum systems the 
quantum body is proven to be convex It is also known, 

that the size of the quantum body may grow with d for two 
parties ([IS [H [12] and for three parties as well 111 . 

In the present paper we wish to further advance the study 
on the shape of the quantum body corresponding to a fixed 
Hilbert space dimension d of bipartite systems. We address 
the problem recently raised by Navascues et al. ITill : Is the 
shape of the quantum body convex for a restricted dimension 
dl Our main result shows that even by four two-outcome 
measurement settings per party the corresponding quantum 
body for a pair of two-dimensional quantum systems (qubits) 
is concave. This result holds for the most general POVM mea- 
surements and for projective measurements as well and can be 
generalized beyond qubits to any dimension d. 

Preliminaries. Let Alice and Bob have two components 
of a compound physical system. Let Alice and Bob choose 
one of a set of rriA and tob two-outcome measurements, re- 
spectively, and let them perform the measurement chosen on 
their respective subsystems. Let us denote the outcome of Al- 
ice's measurement i and Bob's measurement j by Ai = ±1 
{i = 1, . . . , mA) and Bj = ±1 (j = 1, . . . , m^), respec- 
tively. Let us denote the vector having components {Ai), (Bj) 
and (AiBj) for all i and j by x e ]RmA+mB+mAmB^ ^jjere 
(•) denotes the expected value. Vector x may be measured by 
repeating the procedure above on many copies of the system, 
making sure that each pair of measurements is chosen 
to be performed a sufficient number of times to get a satisfac- 
tory statistics. The actual vector one gets will depend on the 
physical system and on the measurement settings the parties 
are allowed to choose from. We note that (Ai) and (Bj) is 
only defined sensibly if the probability of getting a measure- 
ment outcome by one party does not depend on which mea- 
surement the other party has chosen. This is the requirement 
of no-signaling, which is true in both classical and quantum 
physics, and believed to be true in Nature. 

As we have mentioned above, the set of vectors one may 
get when making measurements on systems obeying classical 
physics, or any locally realistic model, is a polytope ll], [^. 
The vertices of the polytope correspond to the deterministic 
situations, when each Ai and Bj has a definite value every 
time it is measured, and the polytope itself is the convex hull 



2 



of these points. A Bell inequality of the form M • x < K 
define an (m^ + + m^jTis — 1) -dimensional hyperplane 
touching the polytope, such that the polytope is on that side of 
the hyperplane, which satisfies the inequality. Tight Bell in- 
equalities are the ones that define the hyperplanes of the facets 
of the polytope. While all vectors x allowed classically may 
be reproduced by measurements on quantum systems, the op- 
posite is not true. In quantum mechanics Bell inequalities may 
be violated, therefore the set of the vectors x allowed is larger 
It is not a polytope, but it is still a convex set . 

In quantum mechanics the components of x may be calcu- 
lated as (Ai) = tv^pAi ® Ib), {Bj) = tr(p/A ® B^), and 
{A^Bj) = iv{pAi ® Bj). Here Ai e Ha (Bj G Hb) is the 
observable corresponding to Alice's (Bob's) measurement i 
(j). Ha (Hb) is the Hilbert space associated to the subsystem 
of Alice (Bob), I a G H-a and Ib Hb are the unity oper- 
ators of the respective Hilbert spaces, and p e Ha <E) Hb 
is the density operator of the physical system. Operators 
Ai = Ia ~ 2P/^ and Bj = Ib - 2P^ have eigenvalues 

±1, while G Ha and P^ £ Hb projection oper- 
ators, whose expected values give the probability of getting 
outcome —1 for the corresponding measurements. If we do 
not confine ourselves to projective measurements, but we al- 
low the more general POVM measurements, then P/^ (Pf) 
will be the POVM element associated to outcome —1 of the 
Alice's (Bob's) measurement i (j), which is not necessarily a 
projector, but any positive operator with eigenvalues between 
and 1. The relation with Ai (Bj) remains the same as above. 

Method. Let us arrange Xqq = 1, and the components 
of X into a matrix of tua + 1 rows and + 1 columns 
as Xm = (A,), Xoj = (Bj) and = {AB^). The 
quantum mechanical expression for these matrix elements is 
Xij = tr {pAi (g) Bj), with the definitions Aq = I a and 
IBq = Ib, here we allowed indices i and j to take value 0. 
Let us restrict the dimensionality of the component Hilbert 
spaces Ha and Hb two. In two dimensions all Hermi- 
tian operators can be written as a real linear combination of 
the three Pauli operators and the unity operator, therefore if 
iriA > 3, there must exist an Ai operator which can be writ- 
ten as a linear combination of the other Ak operators and the 
unity operator with real coefficients. Then the row of the Xij 
matrix depending on that operator can also be written as the 
linear combination of the rows depending on the other Ak op- 
erators and the zeroth row, with the same coefficients. We can 
conclude that correlations described by Xij having more than 
four linearly independent rows can not be reproduced with 
measurements taken on a pair of qubits. We may repeat the 
argument for the columns, too. When we go beyond qubits, 
but still restrict ourselves to finite dimensional Hilbert spaces, 
we can draw a similar conclusion: Xij having more than 
linearly independent rows or columns can not be reproduced 
by measurements performed on systems with no more than d- 
dimensional component Hilbert spaces. This is because any 
d-dimensional Hermitian matrix can be characterized by d 
real (diagonal elements) and d{d — l)/2 complex (nondiago- 
nal elements) numbers, altogether d^ real numbers, therefore, 
no more than d^ of them may be linearly independent with 



real coefficients. The conclusion holds for both projective and 
POVM measurements, as only the hermiticity of the operators 
has been used in the argument. The inclusion of POVM is im- 
portant, because as we will show, if the dimensionality of the 
quantum system is restricted, the quantum body can be larger 
if we allow POVM. 

Main result. Now let us consider the set of x vectors achiev- 
able with measurements on quantum systems of at most d- 
dimensional component Hilbert spaces. The set is not con- 
vex, if there exist points in the vector space that belong to 
the set, but some point on the convex hull of these points 
does not. The latter can be proven by showing that the ma- 
trix Xij corresponding to that point has more than d^ linearly 
independent columns or rows. We will prove below that if 
niA — niB ^ m > d^ and even, the set will not even con- 
tain some vector on the convex hull of points corresponding 
to deterministic cases, that is an element of the local polytope. 
The deterministic cases can obviously be reproduced with any 
physical systems by using degenerate measurements, mea- 
surements with definite outcomes independent of the physical 
system. 

We note that to express the classical or quantum limits on 
results of correlation experiments very often not the vectors 
X, but the vectors p are used, whose components are pA^ , PBj 
and PAiBj , which are the probabilities of getting outcome —1 
for Alice's ith, for Bob's jth, and for both experiments, re- 
spectively. The two approaches are equivalent lUSll . 

Let TUA = mB ~ m he even, and let us take all 
TO!/(m/2)!^ deterministic cases with the Ai being +1 the 
same number of times as it is -1, and Bi — —Ai. Let us call 
the corresponding vectors and matrix elements x'^'^^ and ^^J^ 
((T = 1, . . . , m!/(m/2)!^), respectively. Let x'^^^ and x*^^) be 
defined by = = +1 and A^^'^ = = -1, re- 
spectively. Let us take the following point on the convex hull 
of these vectors: 

x° = !^M^yx('^) + ^fxW+x(-)). (1) 
m m! 2m V / 

(7 

For each deterministic strategy considered in the sum above, 
there is another one with the same weight with all measure- 
ment outcomes having the opposite sign, therefore X°q = 
^oj = = l,...,m). As Xff) = Al'^sf) = -1, 
and = Xfr' = +1' it follows that X° = (2/m) - 1. It 
is easy to see, that if i ^ j, the value of X^^'' = A^'^'^sj'^' = 
-A^'^'^A'f^ is +1 for 2[{m - 2)!]/(?7i/2 - l)!^ cases, and -1 
for the rest of them, and it is obvious that X^/^' = ^1"' = 

J — L 

+1. From these and from Eq. (yj it follows that the non- 
diagonal matrix elements with indices larger than zero are 
Xij = 2/m. The matrix has a nonzero determinant, all 
m + 1 rows and columns are linearly independent, therefore, 
if d^ < m, x° can not be reproduced by measurements per- 
formed on quantum systems with d-dimensional component 
Hilbert spaces. 

Explicit Bell polynomial. Now we will show that all vec- 
tors x°, x'^'^K x(+^ and x'^' considered above belong to a set 
that maximizes a Bell inequality, which can not be violated 
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in quantum mechanics, so they are on the surface of both the 
classical polytope and the quantum set. As the quantum set 
has a multidimensional intersection with the polytope, it fol- 
lows that its surface can not be round everywhere. This fact 
has been also reported recently in the work of Linden et al. 
lfl6ll in the context of distributed computing. The intersection 
has a lower dimensionality than a facet, so the Bell inequal- 
ity is not a tight one. It is a correlation type inequality, that 
is the factors multiplying (Ai) and (Bj) are zero. The Bell 
polynomial is 



m m 



m m 



m , 



i=l 3 = 1 



1=1 3 = 1 



(2) 

where 5ij is the Kronecker delta. To get the maximum value 
of this expression it is enough to consider pure states and pro- 
jective measurements. It is proven in ||l7i1 that for any ob- 
servables A and B in Alice's and Bob's component spaces, 
respectively, and state ifj there exist Euclidean vector a in- 
dependent of 13 and vector h independent of A such that 
{AB) = {'4>\A ® B]^^) = a ■ b. Therefore, we may replace 
{AiBj) with Si ■ bj in Eq. (|2]i, and maximize that expression. 
The vectors ai have to be chosen such that they are parallel 
with the vectors they are multiplied with. Then we get: 



i=l i=l 



(3) 



(4) 



ni ni 



E ■ E ^ ' E 

i=l k=l 



k=l 



We will show that we get the maximum value for B if 



(6) 



k=l 



is true for any i and j. Then one can see from Eqs. ( 1415b that 
li = m/2, and B = This agrees with the upper limit 

this Bell expression may take with quantum measurements, as 
it can be shown analytically making use of semidefinite pro- 
gram ming technique. The actual proof, following Wehner's 
work iflSll is deferred to Appendix [Al 

From Ref. ifTlll it follows, that if itia = tub, and the maxi- 
mum value of the Bell expression can be achieved with all hj 
are linearly independent, than this solution can not be unique. 
The present case is an example for this situation. Equation (|6]l 
has an infinite number of solutions, with hj spanning spaces 
of any dimensionality up to m. An obvious one-dimensional 
solution is when all hj are chosen to be the same unit vector 
h. Then from Eq. Q and U = m/2 it follows that cti = h. 
This arrangement corresponds to the classical deterministic 
strategies of having all measurement outcomes either +1 or 
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FIG. 1 : (Color online) Quantum region in the three-dimensional sec- 
tion spanned by the expectation values (AsBi), {AsB2) and (As) 
as described in the text. The two antipodal apices of the cones, cor- 
responding to and x^, and the equator (in black color) can be 
attained by projective measurements performed on qubits. Whereas, 
any other point on the surface of the cone (represented by red color) 
can be achieved by applying POVM measurements. 



— 1 every time (correlation vectors x*^"*") and x*^")). If m is 
even, there are further one dimensional solutions, with half 
the bj pointing to one direction and the other half pointing to 
the opposite direction. Such a solution corresponds to a de- 
terministic strategy in which Bob has as many measurements 
with a definite outcome of + 1 as ones with an outcome of — 1 , 
and Alice gets the outcome Ai = —Bi for each i (correlation 
vectors x^"^)). From the existence of classical deterministic 
strategies giving the quantum limit for the Bell expression it 
follows that the Bell inequality can not be violated. 

There is an infinite number of solutions of Eq. ^ with bj 
spanning the maximum of m dimensions. An arrangement 
with all hj are orthogonal to each other is one of them. Then 
tti are also orthogonal to each other, and {AiBj) = a; • bj = 
(2/to) — Sij, (see Eq. ([3])). According to Tsirelson's con- 
struction these values can be realized as quantum expecta- 
tion values of ±1 valued observables with a maximally entan- 
gled state of a system of 2™/^ dimensional component Hilbert 
spaces With this state the expectation values (Ai) ~ and 
(Bj) = 0. By choosing the unit vectors corresponding to the 
unity operators in the component Hilbert spaces orthogonal 
to the space spanned by hj, we can get all components of the 
correlation vector as dot products. This correlation vector is 
nothing else than x°, which we have chosen to construct on 
the convex hull of the set of classical deterministic cases x^'^\ 
x^+'' and x^^-* according to Eq. ([T]) as an example that can 
not be achieved with quantum systems of component spaces 
of d < y/m dimensions. Clearly, we could have chosen an 
infinite number of other vectors with the required property. 

POVM versus projective measurements. Now we will show 
that the quantum body with POVM allowed may contain 
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points that can not be achieved with only projective measure- 
ments. Let us consider the quantum body with = 3 and 
niB = 2, restricting ourselves to quantum systems of two di- 
mensional component Hilbert spaces. Let Ai, A2, Bi and 
B2 be the operators, and |?/)) be the pure maximally entangled 
state giving the maximum violation of the CHSH inequality. 
Let the components {Ai), (Bj) and (AiBj) — 1,2) of 
X be derived as the expectation values of the operators above. 
Let the components {A3), {A3B1) and {A3B2) be the expec- 
tation values with ^3 corresponding to a projective measure- 
ment, that is an observable with eigenvalues ±1. Then it can 
be shown that the region allowed for these three components 
are the two antipodal apices of the cones x+ and (when 
A3 = Ia and A3 = —I a, respectively) and the equator of unit 
radius (when A3 has both eigenvalues +1 and —1), shown in 
Fig. [T] To see this, one has to to use the facts that to get the 
other, fixed components of the vector the state must be max- 
imally entangled and the relationship between Bi and B2 is 
also well defined. For example, x'*' = Ax+ + (1 — A)x^ 
(0 < A < 1), a point between the antipodes, can not be 
achieved, as the expectation value of A3 calculated with a 
maximally entangled state can only be ±1 or 0, and when 
it is 0, {A3B1) and {A3B2) can not be at the same time. 
However, we do achieve the point required with the choice of 
A3 = (2A - 1)/^. This operator corresponds to a POVM 
with POVM elements \Ia and (1 — X)Ia associated with the 
+\ and -1 outcome of the measurement, respectively. Sim- 
ilarly, it is easy to prove that all other points within the red 
region shown in Fig.[T]can be attained with POVM. 

Conclusion. We proved that the full set of quantum prob- 
abilities in the bipartite scenario generated either by two- 
outcome projective or by two-outcome POVM measurements 
for any given dimension d is concave. However, one may fur- 



ther ask, whether this fact also holds true for more parties and 
for more than two outcomes. We also proved that the set gen- 
erated by projective measurements may be smaller than the 
one corresponding to the more general POVM measurements. 
In case of two-outcome measurements the maximum viola- 
tion of a Bell inequality with fixed dimensional systems can 
still be achieved with projective measurements fl9Ll20ll . It re- 
mains an open question if this is true in cases of more than 
two outcomes lulll . 

A further question raised by Brunner et al. is that what 
happens, if we restrict ourselves to measurements on a given 
quantum state and look for the set of quantum probabilities 
generated this way. When we limit the dimensionality of the 
Hilbert space, we have shown here that the quantum set is 
concave, by showing that a point on the convex hull of points 
corresponding to deterministic strategies does not belong to 
the set if the number of measurement settings is large enough. 
By restricting ourselves to a particular state, the set can only 
get smaller, while it will still contain the points of determinis- 
tic strategies. 

Finally, it would also be interesting to find out the minimum 
number of settings which generates a quantum body with a 
concave shape for a fixed dimension. In particular, would it 
be possible in the bipartite case to go below four two-outcome 
measurement settings per party by local dimension two in or- 
der to prove concavity of the corresponding quantum body? 
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A. Appendix: Quantum maximum via SDP 

We follow the SDP method put forward by Wehner IHt] 
recently, in order to prove analytically quantum bounds for the 
correlation type Bell inequalities of Eq. (|2]i. Let us consider 
the m X m matrix M with real coefficients 



M^, = 1 - y%, 



(7) 



introduced in Eq. (|2]i. As stated in the main text, the expecta- 
tion values in the polynomial Eq. ^ can be replaced by dot 
product of unit vectors. 



max Mijdi ■ bj 



(8) 



where maximization 

, flm, Ol ..... & 



is taken over all unit vectors 
{fli, . . . , dm, bi, . . . , b,n} € i?^™. As shown by Tsirelson, 
the maximum obtained in this way corresponds to the 
maximum quantum value as well . 

However, the above problem can be formulated as the fol- 
lowing SDP optimization fl^: 



maximize -Tr(rM^) 

subject to r ^ 0, Vz Tu = 1 ■ 



(9) 



Here the matrix W is built up as 



W 



M 
M 



(10) 



and r = ij^ij) is the Gram matrix of the unit vectors 
{ai, . . . , d,m bi, . . . , bm} G i?^™. Denoting the columns of 
the above vectors by V, we can write F = V^V if and only 
if r is positive semidefinite. The constraint Ta = 1, on the 
other hand, owes to the unit length of vectors di and bj . Note, 
that the primal problem defined by ^ is the first step of the 
hierarchy of semidefinite programs given by Navascues et al. 

0B. 

However, one can also define a dual formulation of the SDP 
problem (for an exhaustive review see ll23ll ): 



maximize Tr(diag(A)) 

subject to - + diag(A) >z 0, 



(11) 



where A is a 2m-dimensional vector with real entries and we 
note that this dual problem is just the first step of the hierarchy 
introduced by Doherty et al. ll24ll . 

Let us denote by p* and d* the optimal values for the pri- 
mal and the dual problems, respectively. However, according 
to weak duality, d* > p* fz^ - Thus, in order to prove op- 
timality of the quantum bound one suffices to exhibit a fea- 
sible solution both for the primal (|9]) and for the dual ( fTTT ) 
problem and then show that they are in fact equal to each 
other For this sake let us guess the primal optimum by setting 
(ij, bj = (1, 0, . . . , 0) in dSll with a Bell matrix defined by (Q. 
These vectors correspond to a classical deterministic strategy 
and this solution yields p* = J2Tj ~ iti^/2. 

Similarly, we guess the solution A* = (7n/4)(l, . . . , 1) 
for the dual problem, for which the dual value is d* = 
Tr(diag(A*)) = w? 12. In order to get a feasible solu- 
tion, it remains to check according to (fTTT i whether R = 
— (l/2)VF + diag(A*) >: Ois satisfied. This amounts to prove 
7mi«[-R] > 0, where we use the notation jmin[R] ilmaxiR]) 
for the smallest (largesOeigenvalue of a matrix R. However, 
due to Weyl's theorem ll25ll . for two Hermitian matrices P and 
Q, it holds 7„„„ [P + Q]> jmin [P] + Irnin [Q] ■ For our par- 
ticular case. 



lm.in[R] > Jmin[-^W] + 7„ii„ [diag(A* )] 



m 



(12) 



The eigenvalues of matrix W in the form ( fTOb are given by 
the singular values as ~ VtItJ of matrix M of ^ and their 
negatives. The eigenvalues of M on the other hand are the 
roots of the characteristic polynomial dct(A/ — 7s 1), where 
1 is the m X m unit matrix. In ifTIll we found that the de- 
terminant of an m X m matrix with diagonal elements p and 
non-diagonal elements g is [p + (m — l)q](p — q)"''^^. By 
inserting p ^ \ — m/2 — and q = 1 into the determinant 
above, we obtain the roots 7^ ~ ±m/2. This result implies 
7maa:[W^] = ^. By Substituting this value into ( fT2l i we get 
Imin [P\ > 0. This implies that this solution for d* is feasi- 
ble, and recalling the guessed solution p*, we have d* ~ p* ■ 
Thus the maximum quantum value of the Bell polynomial M 
defined by Eq. ([2]) is equal to m? /2, which can be achieved by 
classical means as well. 



